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ABSTRACT 

We study the distribution of projected ellipticity n(e) for galaxies in a sample of 
20 rich (Richness ^ 2) nearby (z < 0.1) clusters of galaxies. We find no evidence of 
differences in n(e), although the nearest cluster in the sample (the Coma Cluster) is 
the largest outlier (P(same) < 0.05). We then study n(e) within the clusters, and find 
that e increases with projected cluster-centric radius R (hereafter the e-R relation). 
This trend is preserved at fixed magnitude, showing that this relation exists over and 
above the trend of more luminous galaxies to be both rounder and more common 
in the centres of clusters. The e-R relation is particularly strong in the subsample 
of intrinsically flattened galaxies (e > 0.4), therefore it is not a consequence of the 
increasing fraction of round slow rotator galaxies near cluster centers. Furthermore, 
the e-R relation persists for just smooth flattened galaxies and for galaxies with de 
Vaucouleurs-like light profiles, suggesting that the variation of the spiral fraction with 
radius is not the underlying cause of the trend. We interpret our findings in light of 
the classification of early type galaxies (ETCs) as fast and slow rotators. We conclude 
that the observed trend of decreasing e towards the centres of clusters is evidence for 
physical effects in clusters causing fast rotator ETGs to have a lower average intrinsic 
ellipticity near the centres of rich clusters. 
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1 INTRODUCTION 

Early Type Galaxies (ETGs) account for half of the 
stellar mass in the local Universe (Renzini 2006). They are 
traditionally divided in two subclasses: Elliptical (E) and 
lenticular (SO) galaxies. To first order, Es have smooth, 
single component light profiles, while SOs present both 
a central bulge and an extended stellar disc. ETGs are 
more common in clusters of galaxies, and the morphology- 
density relation (T-E, Dressier 1980) illustrates the 
effect of the local environment on the formation and 
evolution of these galaxies (Whitmore et al. 1993, argue 
however that the morphology correlates better with the 
cluster-centric radius). Dressier et al. (1997) used HST 
photometry to show that cluster Es are already in place 
at redshifts z ~ 0.5, while the fraction of SOs is lower 
than in the local Universe. Therefore Es form earlier 
than SOs, which arise from infalling Late Type Galaxies 
(LTGs) mostly between 2 « 0.5 and 2 = 0 (e.g. Vulcani 
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et al. 2011, but see Holden et al. (2009) for a different view). 

The division between Es and SOs presents however 
a number of problems. Observationally it is difficult to 
distinguish Es from close-to-face-on SOs, and morphological 
catalogues might be biased in this sense (van den Bergh 
1990). Galaxies classified as Es are often found to contain 
disc components (Kormendy & Djorgovski 1989). Rix 
& White (1990) further demonstrated how discs can go 
undetected even in local galaxies. All but the brightest Es 
and SOs in the Coma Cluster form a family with continuous 
bulge-to-disc light ratios (B/D, Jprgensen and Franx 1994). 
More recently, the SAUR0N survey (de Zeeuw et al. 2002) 
used Integral Field Spectroscopy (IFS) to investigate the 
stellar kinematics of ETGs. They identified two dynamical 
classes within ETGs: Fast Rotators (FRs) and Slow Rota¬ 
tors (SRs). The former exhibit large scale rotation patterns 
typical of a disc origin, while the latter have little to no 
rotation (Emsellem et al. 2007; Cappellari et al. 2007). 
Importantly both SRs and FRs are found amongst both Es 
and SOs. The volume-limited ATLAS 3D survey (Cappellari 
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et al. 2011a) determined that 66% of the local Es are FRs 
(Emsellem et al. 2011), and a new classification paradigm 
has been invoked (Cappellari et al. 2011b; Kormendy & 
Bender 2012). 

In particular, FR ETGs form a sequence of increasing 
disc fraction, parallel to that of spiral galaxies (the comb 
diagram, Cappellari et al. 2011b). The emerging picture is 
that - kinematically - FR ETGs are much more similar to 
spiral galaxies than they are to SR ETGs. Unfortunately 
IFS observations, that are necessary to tell apart FRs from 
SRs, are time consuming when compared to photometry. 
Even the largest upcoming IFS surveys (SAMI, MaNGA; 
Groom et al. 2012; Bundy et al. 2015) have sample sizes 
which are several orders of magnitude smaller than current 
state-of-the-art photometric surveys, like SDSS. Here we 
investigate whether we can take advantage of currently 
available photometric samples to investigate the proper¬ 
ties of ETGs in the framework of the FRs/SRs classification. 

Using photometry alone, ETGs can be characterised by 
their projected ellipticity e = 1 — q, where q is the appar¬ 
ent axis ratio. The observed value of e depends both on the 
distribution of the orbits (which determines the intrinsic el¬ 
lipticity) and on the inclination of the galaxy on the plane 
of the sky. In general it is not possible to infer the 3-D struc¬ 
ture of an individual galaxy from its photometry alone, but 
we can study intrinsic shapes statistically (Sandage et al. 
1970; Lambas et al. 1992). These studies showed that the 
distribution of ellipticity n(e) for Es and SOs is different, 
with Es on average rounder than SOs. More recently, Wei- 
jmans et al. (2014) used the ATLAS 30 sample of ETGs to 
show that the intrinsic shape of FRs (both Es and SOs) is 
consistent with that of spiral galaxies, with an average axis 
ratio q = 0.25±0.01. This reinforces the view that FR ETGs 
and spiral galaxies form a family of intrinsically flat stellar 
systems. 

Kuehn & R.yden (2005, KR05) studied the relation 
between q and the local environment using a magnitude 
limited sample (r < 17.77 mag). They found that galaxies 
with different magnitude and light profile exhibit different 
trends of t with the local number density of galaxies. 
Galaxies characterised by a de Vaucouleurs light profile 
(de Vaucouleurs 1948) are rounder in denser environments, 
regardless of their luminosity. Galaxies with exponential 
profiles show two opposite tendencies, based on their 
absolute magnitude. Galaxies fainter than M r = —20 mag 
are rounder in denser environments, while more luminous 
galaxies tend to become flatter at higher density. These 
results suggest that the environment can affect the shape 
of galaxies. Intriguingly, the trend observed by KR05 is 
strongest when the density is measured inside an aperture 
with diameter 2 /r~ 1 Mpc, which roughly corresponds to the 
size of a cluster of galaxies. 

Here we propose to investigate n(e) for cluster galaxies, 
in the framework of the SR/FR paradigm. In particular, 
since there are no SRs with e Js 0.4 (Emsellem et al. 2011), 
we can remove them from any photometric sample. In the 
next section we introduce our sample of local cluster galax¬ 
ies, and proceed to show the results of the analysis. We then 
discuss possible sources of bias, and whether the observed 


relation is a consequence of previously known trends. We 
conclude with a summary of our results. 


2 DATA AND SAMPLE 

To assess the effect of the environment on the apparent el¬ 
lipticity of galaxies we study a sample of rich nearby clusters 
of galaxies with SDSS DR10 data (Ahn et al. 2014). At low 
redshift (z < 0.04), the most complete catalogue of clusters 
is that of Abell et al. (1989), from which: 

• We discarded all the clusters with 2 > 0.1, or with no 
redshift measurement available. 

• We discarded all clusters with richness 1Z < 2 because 
we expect any signature of the cluster environment to be 
strongest in rich systems. 

This left 20 clusters (see Table 1). For each of them, we 
retrieved SDSS photometry as follows. 

• We determined the centre of the cluster as the coor¬ 
dinates of the brightest galaxy. When the second brightest 
galaxy falls within 0.5 mag from the brightest, we used the 
midpoint between the brightest and second brightest as the 
centre of the cluster (e.g. Abell 1656, Abell 1367). 

• We queried the table “Galaxy” from the SDSS DR10 
database, retrieving all the galaxies within a projected ra¬ 
dius of 1.5 Mpc from the cluster centre. The angular dis¬ 
tance was derived by assuming Planck Cosmology (fl m = 
0.32,11a = 0.68, ho = 0.67, Ade et al. 2014). We assume 
zero peculiar velocity for all the clusters. 

• We removed all galaxies with a radius deVRad 1 ^ 0.4” 
(corresponding to 0.5 kpc at 2 = 0.1). This constraint filters 
out artefacts and measurements with bad photometry, which 
are otherwise still present in the table Galaxy. 

• We further removed objects that are classified as stars 
in the g'- and i'-band (SDSS uses only U-band for this classi¬ 
fication). This condition filters out 31 objects. Objects with 
negative or zero errors in the photometry were also elimi¬ 
nated (21 objects in total). 

• We applied a cut in absolute magnitude at M r = 
— 18.0 mag. This is well above the r'-band completeness 
limit of SDSS at redshift 2 = 0.1 (the maximum distance 
modulus for our sample of clusters is « 38 mag). A more gen¬ 
erous magnitude cut increases the scatter in the RS, which 
is undesirable. We do not apply a k-correction, but this does 
not affect our results. 

For each galaxy we retrieved the SDSS r'-band Model- 
Mag as a measure of magnitude and g' — r' colours using 
SDSS AperMag. 

Ancillary data consist of SDSS DR10 spectroscopic red- 
shifts for a sample of galaxies in the cluster Abell 1656 (the 
Coma Cluster), and of Galaxy Zoo 2 (GZ2, Willett et al. 
2013) morphological classifications, when available. 


1 Column names from the SDSS table are reported in italics. 
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Figure 1. Comparison between e mo del and e m oments> measured in 
the r'-band as described in the text. The solid red line is the best- 
fit linear relation, while the black dashed line is the 1:1 relation. 
Emodel is on average higher than e m0 ments. 


2.1 Projected ellipticity 

SDSS offers two different measures of ellipticity. The model 
ellipticity is defined as 1 — b/a, where b/a is the apparent 
axis ratio of the r'-band best fit model (i.e. b/a= deVAB_r 
if fracDeV-r ^ 0.5 and b/a = expAB_r if fracDeV-r< 0.5). 
The SDSS pipeline automatically corrects the model axis 
ratio for the effect of the Point Spread Function (PSF; see 
Section 4.1 for a discussion). 

An alternative measure of the ellipticity is to use the 
adaptive second order moments of the surface brightness 
/. This method uses a Gaussian weight function adaptively 
matched to the size and shape of the galaxy being measured 
(Bernstein & Jarvis 2002; Hirata & Seljak 2003). 

Comparing the model ellipticity to the mo¬ 
ments ellipticity in the r'-band, we find that in gen¬ 
eral the former is larger. The best linear fit yields 
emodei = (1-039 ± 0.004)e mome „ts + (0.025 ± 0.001) (Fig. 1). 
As in general e increases with radius within galaxies, 
fmodei is better suited for the outer region of galaxies. 
This is reinforced by the fact that most of the cases where 
^moments Cmodei are close-to-face-on barred galaxies. 
We are primarily interested in the effect of the cluster 
environment on the shape of galaxies. Since this effect is 
larger at larger galactic radius, we choose to use e mo dei as 
a measure of the shape of galaxies. See Appedix A for a 
comparison with the results using e momen ts. 

We investigate the variation of e as a function of two 
different tracers of the environment density, the projected 
cluster-centric radius R and the projected number density of 
galaxies E 3 . The latter was measured for each galaxy inside 


the circle on the sky comprising its three closest neighbours 
(e.g. Cappellari et al. 2011b). 

2.2 Red Sequence determination 

For each cluster we constructed the g' — r' vs r' Colour- 
Magnitude Diagram (CMD). In order to reject any outliers 
and artefacts we eliminated all the points in the CMD with 
g' — r' (/ [0.0, 2.5] and galaxies with errors in g — r' greater 
than 0 . 1 . 

We identified the Red Sequence (RS) using the Gaus¬ 
sian mixture model of Houghton et al. (2012). The algorithm 
fits two superimposed Gaussian models: one for the RS and 
one for the underlying galaxy distribution, including back¬ 
ground and foreground objects. The mean (g 1 — r') of the RS 
Gaussian is allowed to vary linearly with the magnitude r', 
i.e. (g 1 —r')=m (r' — 16) + c, where m and c are constants 
to be determined. The amplitude and the dispersion A of 
the Gaussian are held constant with r'. For the background 
Gaussian, all parameters are independent of r'. 

The algorithm determines the most likely values of m, 
c and A. In addition, it returns the probability P(RS) that 
each galaxy on the CMD belongs to the RS, and we define 
galaxies with P(RS) > 0.05 to be members of the RS. This 
corresponds to a 2<r selection. After applying this procedure 
to all the 20 clusters in the sample, the average results are: 
(m) = -0.037 ± 0.006 and (A) = 0.035 ± 0.007. 

In Fig. 2 we show the CMD of two clusters, Abell 
1650 at 2 = 0.08 (top panel) and Abell 2199 at z = 0.03 
(bottom panel). Each black dot represents a galaxy from 
the SDSS catalogue, but only the dots with an overlaid 
circle have been considered for the purpose of the RS 
determination. Red circles mark RS galaxies, while green 
circles mark non-RS galaxies, as determined by their value 
of P(RS). The vertical dashed line is the magnitude cut for 
the cluster, while the solid green line represents the best fit 
to the RS, i.e. the line of equation g' — r' = m (r' — 16) + c. 
In the top left box of each panel we report the best fit value 
of m, c and A with the corresponding 1 a confidence interval. 

We constructed a Reference Cluster (RC) as the union 
of all the clusters in the sample, and further constructed a 
reference RS as the union of all the RS’s of all the clusters. 
The RC consists of 9052 galaxies, the reference RS sample 
counts 5175 entries (see last row in Table 1, columns 7 and 
8). 

2.3 Redshift selection 

A colour selection can be useful to reject interlopers, but it is 
not as reliable as a spectroscopic selection. Individual cluster 
galaxies can be significantly bluer than allowed by the RS 
(for instance, if they are undergoing star formation). On the 
other hand, a number of physical processes can affect the 
colour of interlopers, placing them close or on the RS (dust 
reddening, star formation). In order to assess the number 
of interlopers still present in the RS sample, we used spec¬ 
troscopic data from SDSS DR10. Unfortunately only Abell 
1656 had enough candidate galaxies with spectra, therefore 
the spectroscopic sample is limited to one cluster. We re¬ 
trieved the redshift a for all the galaxies in the photomet¬ 
ric sample, and rejected all the galaxies with 2 < 0.01 or 
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Figure 2. Two example CMDs for Abell 1650 at 2 = 0.08 (top) 
and Abell 2199 at z = 0.03 (bottom). Galaxies are marked by 
black dots. Red and Green circles are overlaid on RS galaxies 
and to non-RS galaxies. Naked dots denote galaxies with large 
errors that were not used in the fit. The best fit to the RS is the 
green line, while the region shaded in red represents the intrinsic 
scatter. The vertical dashed line is the adopted magnitude cut. 
The best fit parameters are reported in the top left box (m is the 
slope, c the intercept and A the intrinsic scatter). 

z > 0.04 (Price et al. 2011). The sample for Abell 1656 is 
reduced from 468 to 387 galaxies, while the RS is reduced 
from 362 to 333 members. The number of interlopers is 61 
and 21 for the full sample and for the RS respectively, cor¬ 
responding to a fraction of 0.14 and 0.06. We then assumed 
that all of the galaxies with no redshift measurement in the 
database are interlopers too (20 and 8 galaxies for the full 
sample and for the RS respectively). This gives a maximum 
fraction of interlopers of 0.18 and 0.08. 


3 RESULTS 

We start this section studying if n(e) varies between different 
clusters, and the properties of the RS sample (Section 3.1). 
We then study the redshift selected sample of Abell 1656 to 
study the effect of interlopers (Section 3.2). In Section 3.3 
we consider again the full sample, where we repeat the anal¬ 
ysis for intrinsically flat galaxies. We then study the relation 
between e and r'-band luminosity (Section 3.4) and conclude 
this section by studying a subsample of morphologically se¬ 
lected galaxies (Section 3.5) and splitting our sample ac¬ 
cording to the shape of the luminosity profile (Section 3.6). 



e 

Figure 3. Distribution of projected e for the Reference Cluster 
(solid green line) and the reference RS (solid red line). The dashed 
lines show the distribution for Abell 1656, with the same colour 
coding. According to a KS test (Table 1), Abell 1656 has a dif¬ 
ferent e distribution than the Reference Cluster. This is true also 
for the RS. The distributions are normalised so that the integral 
is unity, but the values of n for Abell 1656 are plotted 0.3 higher 
for display purposes. 


3.1 Full sample 

In Fig. 3 we plot n(e) of the RC (solid green line), alongside 
the distribution of the reference RS (solid red line). Each 
distribution has been normalised so that its integral between 
e = 0 and e = 1 is unity. For both the RC and the RS we 
find an average ellipticity (e) = 0.38 ± 0.01 and a standard 
deviation <r e = 0.21 ± 0.01. 

We then used a Kolmogorov-Smirnoff (KS) test on the 
null hypothesis that any individual cluster in the sample 
has the same n(e) as the RC, and that any of the cluster 
RS has the same n(t) as the reference RS. The results 
are summarised in Table 1. For each cluster we list the 
number of galaxies found inside the adopted SDSS aperture 
(column 7) and the number of galaxies on the RS (column 
8). We then report the results of the KS test: for each 
cluster we list the probability that its n(e) is the same as 
the RC (column 11), and that n(t) of its RS is the same 
as the reference RS (column 12). We have only one entry 
with P < 0.05 for both column 11 and 12 (Abell 1656 and 
Abell 16 respectively). If we assume that n(e) is the same 
between all clusters, we expect the results of the KS test to 
be distributed uniformly between 0 and 1. We performed 
an Anderson-Darling (AD, Anderson & Darling 1954) 
uniformity test for the null hypothesis that the values of 
columns 11 and 12 of Table 1 are drawn from the uniform 
distribution. The AD values are 0.48 and 0.41 for columns 
11 and 12 respectively, corresponding to p-values 0.25 
(Rahman et al. 2006). Therefore there is no evidence of any 
difference between the n(e) of our sample of clusters. 

Next we studied the variation of t with the cluster¬ 
centric radius R. 

In Fig. 4 we plot e for nine bins in log R. The green di¬ 
amonds mark each bin from the parent sample in the CMD, 
and the red circles mark e for the galaxies on the RS. At 
the top of each panel, the numbers in parentheses are the 
numbers of RS galaxies in each bin (we eliminated any bin 
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Figure 4. The projected ellipticity e of galaxies as a function 
of the projected cluster-centric radius. The Reference Cluster is 
shown at the top, followed by two example clusters: Abell 1650 
(centre) and Abell 2142 (bottom). The green diamonds represent 
all galaxies in the cluster; red circles represent RS galaxies only. 
The boxes in the bottom left corners report the best fit slope to 
the points and associated errors. The reference sample, as well 
as Abell 1650, show a trend of e with log R, which we call the 
e-R relation. In Abell 2142 the trend is not significant (the slope 
is less than three standard deviations above 0). The numbers in 
parentheses are the number of galaxies in each bin of log R for 
the RS. 

with less than three galaxies). The errorbars were derived 
assuming that for each bin in log R, the distribution n(e) 
is flat between e = 0 and e = 1, and zero otherwise. Such 
a distribution has standard deviation l/%/l2, therefore the 
error on the mean e for each bin is 1/\/l2 N, where N is the 
number of galaxies in that bin. Since n{e) is hardly flat, and 
does not extend to e = 1 (Fig. 3), it follows that our estimate 
of the error is a conservative one. The box in the bottom left 
corner reports the slope (with errors) of the best linear fit, 
using least squares minimization. The uncertainties on the 
slopes were derived from the errorbars, therefore they too 
represent a conservative estimate. The top panel depicts the 
results for the RC: both the parent and RS sample show a 


clear trend of e with log R, in that the best-fit slope is more 
than three standard deviations away from 0. For brevity, we 
call the observed trend the ellipticity-radius relation (e-R). 
The central panel shows the same plot for Abell 1650, a clus¬ 
ter where the e-R trend is clearly present. In contrast, Abell 
2142 (bottom panel) is an example of a cluster where the re¬ 
lation is not detected, i.e. the measured slope is within three 
standard deviations of 0. For the parent samples, four clus¬ 
ters have a slope that is more than three a from 0, including 
Abell 1650. Nine clusters have slopes between one and three 
a above 0, and seven clusters (including Abell 2142) have a 
slope that is within one a from 0 (this includes Abell 16 and 
Abell 168, which have negative slopes). The breakdown is 
4/9/7. For the RS samples the breakdown is 5/10/5 (Abell 
16 and Abell 168 still have negative slopes). In all clusters 
the best-fit slope of the parent sample and of the RS sample 
are statistically consistent to the level of three a. 

To estimate the probability that our result arises from 
chance, we assume Gaussian errors on the best-fit slope and 
use the binomial distribution: 

f{k;n,p)= ^jp k (l-p) n ~ k (1) 

The probability that out of 20 best-fit slopes 13 are more 
than one sigma above 0 is /(13, 20,0.16) fs 1 x 10 -6 . We 
therefore conclude that the observed e-R relation is a real 
effect. 

We find no correlation between the value of the best fit 
slope and global cluster parameters, like X-ray temperature 
or X-ray luminosity (Ebeling et al. 1996), Richness (Abell 
et al. 1989) or Bautz-Morgan type (Bautz & Morgan 1970). 

3.2 Spectroscopic sample for Abell 1656 

In Section 2.3 we presented a sample of redshift selected 
cluster member for Abell 1656. The e-R relation stays the 
same for the original sample (Fig. 5, top panel) and for the 
redshift selected sample (Fig. 5, bottom panel). The signifi¬ 
cance of the trend increases for both the parent sample and 
the RS, but the trends are statistically consistent. If any¬ 
thing, the slope of the best-fit relation increases going from 
the original sample to the spectroscopic sample. 

3.3 Flat galaxies 

We repeated the analysis for the subset of flat galaxies. Em- 
sellem et al. (2011) showed that there are no SRs flatter than 
e > 0.4 (we treat double sigma galaxies as FRs, Krajnovic 
et al. 2011). Therefore we isolate a sample of intrinsically 
flat galaxies by selecting e ^ 0.4 2 . We performed a KS test 
comparing the flat galaxies in the RC to those in each sin¬ 
gle cluster. We then repeated the test for RS galaxies only. 
The results are listed in Table 1 (columns 13 and 14). We 
find that Abell 1656 fails both tests and Abell 1650 fails one 
test. To test the null hypothesis that the values in columns 

2 Notice that this biases our sample against FRs with lower in¬ 
trinsic ellipticity, because the range of possible inclinations that 
satisifies the constraint e 0.4 decreases with decreasing intrinsic 
ellipticity. 
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Figure 5. Same as Fig. 4, for the Abell 1656 sample (top panel) 
and for the sample of redshift selected cluster members (bottom 
panel; see Section 2.3). The trend of e with cluster centric radius 
R observed in Abell 1656 is not an effect of interlopers: the slope 
of the best-fit linear relations are statistically consistent between 
the original sample and the redshift selected samples (boxes in 
the bottom left corner of each panel). 


13 and 14 are distributed uniformly between 0 and 1, we use 
again the AD uniformity test. The resulting values are 0.53 
and 1.84 for columns 13 and 14 respectively, corresponding 
to p-values 0.25 and 0.10 < p — value < 0.15. We in¬ 
fer that there is no conclusive evidence to suggest that the 
distribution of e is different between different clusters, as 
we cannot reject the null hypothesis with confidence greater 
than P — 0.01. 

However, by looking at Abell 1656 in Table 1, we see 
that this cluster scores very low confidence in three out of 
four tests (columns 11, 13 and 14). Since these values are 
not independent, as a conservative estimate we adopt the 
largest one as the probability that the n(e) of Abell 1656 is 
the same as the n(e) of the RC, P(same) ^ 0.05. 

We now look into the relation of e with the cluster cen¬ 
tric radius R. The top panel of Fig. 6 shows the t-R relation 
for the subsample of flat galaxies. The magenta circles are 
galaxies with t < 0.4 (binned in log R), and the blue el¬ 
lipses are galaxies with t ^ 0.4. The trend is observed for 
both the round and flat subsamples, however the t-R re¬ 
lation for flat galaxies (0.035 ± 0.005) is both steeper and 
and statistically more significant than that of round galaxies 
(0.014 ± 0.003). The bottom panel of Fig. 6 shows the same 
relation for galaxies on the RS: the best fit slope for flat 
galaxies is consistent with the previous one. The change in 



log R [kpc] 

Figure 6. The projected ellipticity e of galaxies as a function 
of the projected cluster-centric radius. The Reference Cluster is 
shown at the top, the reference RS sample in the bottom panel. 
In each plot, the blue ellipses represent galaxies that appear flat 
on the sky (e js 0.4). The magenta circles represent galaxies that 
appear round on the sky (e < 0.4). The blue and magenta lines are 
least squares best fits to the data; their slopes are reported in the 
bottom left corner of each figure. The number in parentheses are 
the number of galaxies in each bin of log R for the flat subsample. 

mean t for flat galaxies is A t~ 0.1 over the observed range 
in log R. 

If we repeat our analysis substituting log R with log E 3 , 
we find equivalent results: e decreases with log E 3 . This is 
a consequence of the anticorrelation between R and E 3 . At 
this stage it is not possible to disentangle which relation is 
more fundamental. 

3.4 Dependence on luminosity 

In Fig. 7 we show how the RC populates the t vs M r -space. 
The colourbar represents the log number density of galax¬ 
ies. We notice that there is an excess of round, luminous 
galaxies at M r « —23 mag (Fig. 7, bottom left corner). 
In contrast, there are only 31/9052 galaxies with e ^ 0.4 
and M r ^ —22 mag, and only 2/9052 galaxies with t ^ 0.4 
and M r ^ —23 mag. It is possible that a number of these 
are foreground interlopers. The brightest of these galaxies 
is however the cD galaxy of Abell 1650, which is genuinely 
flat. For M r < —22 mag, the contour lines indicate a gradual 
increase of e with decreasing r'-band luminosity. 

In Fig. 8 we show t for nine bins in M r , for the RC as 
well as for the reference RS (top panel). For galaxies with 
—22 mag < M r < —19 mag, there seems to be no depen¬ 
dence of t with M r . However, galaxies more luminous than 
~ —23 mag are on average rounder. In the middle and bot¬ 
tom panel of Fig. 8 we plot t against M r for the RC and for 
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Table 1. The local sample of clusters and results of the KS tests. 


ACO 

RA 

DEC 

z 

r 

R 

N parent 

Nrs 

N^o.4 

N RS,e^0.4 

P parent 

Prs 

Pe^0.4 

PRS,e^0.4 


(J2000) 

(J2000) 


(mag) 

(arcsec) 









(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(ii) 

(12) 

(13) 

(14) 

16 

0:16:46.30 

6:44:39.84 

0.084 

19.88 

16.09 

416 

185 

187 

99 

0.875 

0.031 

0.752 

0.656 

168 

1:15:9.79 

0:14:50.64 

0.045 

18.48 

28.52 

330 

181 

148 

84 

0.959 

0.846 

0.854 

0.555 

1035 

10:32:7.20 

40:12:33.12 

0.080 

19.77 

16.80 

461 

194 

207 

94 

0.574 

0.459 

0.147 

0.212 

1186 

11:13:51.36 

75:23:39.84 

0.079 

19.75 

16.95 

311 

152 

128 

63 

0.395 

0.950 

0.533 

0.408 

1190 

11:11:46.32 

40:50:41.28 

0.079 

19.75 

16.89 

358 

224 

162 

112 

0.754 

0.159 

0.365 

0.211 

1367 

11:44:29.52 

19:50:20.40 

0.021 

16.83 

58.29 

284 

206 

118 

87 

0.161 

0.219 

0.149 

0.076 

1650 

12:58:46.32 

-1:45:10.80 

0.085 

19.90 

15.97 

461 

290 

220 

138 

0.239 

0.649 

0.006 

0.117 

1656 

12:59:48.72 

27:58:50.52 

0.023 

16.99 

54.13 

468 

362 

197 

153 

0.046 

0.105 

0.004 

0.010 

1775 

13:41:55.68 

26:21:53.28 

0.072 

19.54 

18.38 

456 

234 

202 

102 

0.824 

0.962 

0.714 

0.692 

1795 

13:49:0.48 

26:35:6.72 

0.062 

19.20 

21.14 

421 

220 

176 

85 

0.240 

0.340 

0.307 

0.633 

1904 

14:22:7.92 

48:33:22.32 

0.071 

19.49 

18.76 

442 

252 

204 

118 

0.902 

0.314 

0.912 

0.233 

2029 

15:10:58.80 

5:45:42.12 

0.077 

19.68 

17.43 

558 

285 

260 

129 

0.381 

0.825 

0.535 

0.578 

2065 

15:22:42.72 

27:43:21.36 

0.072 

19.54 

18.40 

567 

305 

263 

140 

0.414 

0.483 

0.206 

0.424 

2142 

15:58:16.08 

27:13:28.56 

0.090 

20.04 

15.10 

646 

405 

281 

158 

0.708 

0.092 

0.442 

0.426 

2151 

16:5:14.88 

17:44:54.60 

0.037 

18.04 

34.41 

355 

189 

168 

86 

0.305 

0.460 

0.254 

0.143 

2199 

16:28:36.96 

39:31:27.48 

0.030 

17.59 

41.80 

327 

201 

145 

80 

0.847 

0.527 

0.827 

0.529 

2244 

17:2:43.92 

34:2:48.48 

0.099 

20.27 

13.82 

556 

301 

247 

133 

0.573 

0.873 

0.172 

0.361 

2255 

17:12:30.96 

64:5:33.36 

0.081 

19.80 

16.61 

636 

398 

294 

186 

0.346 

0.262 

0.727 

0.797 

2256 

17:3:43.44 

78:43:2.63 

0.060 

19.12 

21.82 

528 

341 

217 

138 

0.317 

0.137 

0.752 

0.412 

2670 

23:54:10.08 

-10:24:18.00 

0.076 

19.66 

17.56 

471 

250 

206 

109 

0.474 

0.676 

0.974 

0.169 

Reference 

- 

- 

- 

- 

- 

9052 

5175 

4030 

2294 

- 

- 

- 

- 


Column (1): cluster ID (Abell et al. 1989). Column (2): right ascension in degrees and decimal. Column (3): declination in degrees and 
decimal. Column (4): cluster redshift from Ebeling et al. (1996), or from Abell et al. (1989) when available. Column (5): magnitude cut 
adopted, SDSS r-band. Column (6): cut in the projected distance, corresponding to 1.5 Mpc. Column (7): number of galaxies in the 
parent sample. Column (8): number of galaxies in the RS. Column (9): number of flat (e ^ 0.4) galaxies in the parent sample. Column 
(10): number of flat (e ^ 0.4) galaxies in the RS. Column (11): probability that the parent sample has the same e distribution as the 
corresponding Reference Cluster. Column (12): probability that the RS has the same e distribution as the corresponding Reference 
Cluster. Column (13): probability that the sample of flat (e ^ 0.4) galaxies has the same e distribution as the corresponding Reference 
Cluster. Column (14): probability that the sample of flat (e ^ 0.4) RS galaxies has the same e distribution as the corresponding 
Reference Cluster. Probability values P < 0.05 are highlighted in boldface characters (see Section 3.1 and Section 3.3). 



H [mag] 

Figure 7. The distribution of the galaxies (Reference Cluster) in 
the magnitude-e plane, n is the number density of galaxies in the 
e — M r space. There is an excess of round, luminous galaxies (M r 
< —22 mag, bottom left corner). 


the reference RS respectively, but this time we divide each 
sample in two two subsets at e = 0.4. Flat galaxies do not 


Table 2. The e-R relation at fixed r'-band absolute magnitude 
(see Section 4.2). 


M r 

(1) 

Ulparent 

(2) 

UT-parent 

^parent 

(3) 

“RS 

(4) 

m rs 

&RS 

(5) 

Mr -22 

0.094 

2.3 

0.089 

2.6 

-22 <Mr4 -21 

0.024 

3.5 

0.028 

2.2 

-21 < Mr ^ -20 

0.038 

1.5 

0.058 

2.1 

-20 < M r ^ -19 

0.034 

2.9 

0.060 

4.4 

-19 < Mr s; -18 

0.067 

6.5 

0.090 

10.0 


Column (1): magnitude interval (r'-band absolute magnitude). 
Column (2): the measured slope for the e-R relation. Column (3): 
significance on the measured slope for the e-R relation, (in units 
of the uncertainty cr). Column (4): the measured slope for the 
e-R relation for RS galaxies only. Column (5): significance on the 
measured slope for the e-R relation for RS galaxies only (in units 
of the uncertainty cr). 

present a trend of e with magnitude, and the slope of the 
best-fit relation is statistically consistent with zero. 

Next we studied the e-R relation by dividing the sample 
in luminosity bins. The e-R relation persists, but decreases 
in significance when we decrease the sample size. The results 
are listed in Table 2. We report the adopted magnitude in¬ 
terval (column 1), the measured best fit slope (column 2) 
and the significance of the e-R relation in that interval (col- 
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M r [mag] 

Figure 8. Top panel: e vs M r for the Reference Cluster (green 
diamonds) and for the RS (red circles). Galaxies more luminous 
than M r RS —23 mag are on average rounder than the galaxies 
with M r > —23 mag. Middle and bottom panel: same as Fig. 6, 
but showing e as a function of M r for e < 0.4 and for e ^ 0.4 . 
When we consider the subsample of flat galaxies (e > 0.4), there 
is no trend of e with M r . This is quantified by the slope of the 
best-fit linear relation, which is statistically consistent with zero 
(bottom left corner of each panel. 

umn 3). We list also the best fit slope for the RS (column 
4) and the significance of the t-R relation for the RS (col¬ 
umn 3). The significance of the t-R relation is the ratio of 
the best fit slope m by its uncertainty a. In general, the t- 
R relation is recovered in each magnitude bin, however the 
strongest slope and the highest significance are found for the 
least luminous galaxies. 

3.5 Galaxy Zoo 2 Morphological Selection 

GZ2 is a citizen science project providing morphologi¬ 
cal classifications for rs 300, 000 galaxies drawn from the 
SDSS sample. Cross-correlating their catalogue with our 
sample we find 1425 matches (rs 16%). We used debi- 
ased likelihoods to identify smooth objects (Psmooth = 



log R [kpc] 

Figure 9. Same as Fig. 6, but for the subsample of cluster galax¬ 
ies morphologically classified as smooth in Galaxy Zoo 2. In each 
plot, the blue ellipses represent galaxies that appear flat on the 
sky (e > 0.4). The magenta circles represent galaxies that ap¬ 
pear round on the sky (e < 0.4). The blue and magenta lines are 
least squares best fits to the data; their slopes are reported in 
the bottom left corner of each figure. The number in parenthe¬ 
ses are the number of galaxies in each bin of log R for the flat 
subsample. The top panel shows the results for the subsample of 
galaxies with P smoo th > 0.8, but it does not have enough galaxies 
to constrain the slope. The bottom panel shows the results whith 
a more generous selection P smoo th >0.5 and P ap iral <0.5. 

t01-smooth-or_features-a01-smooth_debiased). Objects with 
P B mooth > 0.8 correspond to ETGs in the Hubble classifi¬ 
cation (Willett et al. 2013). However, this selection leaves 
us with 252 galaxies, of which only 50 are flat (e ^ 0.4). 
This is insufficient to constrain the slope of the t-R relation 
(m = 0.065 ± 0.036, see Fig. 9, top panel). 

If we relax the above condition and use instead 
Psmooth > 0.5 and P sp irai < 0.5, we find 676 galaxies (and 
204 flat galaxies). In the bottom panel of Fig. 9 we show 
that the t-R relation still holds for flat, smooth galaxies, 
albeit only to the 2 a level (the best-fit linear slope is 
m = 0.033 ± 0.016). This is statistically consistent with the 
trend without morphological selection (m = 0.035 ± 0.006). 
Therefore it appears that the spiral contamination is not the 
driving mechanism behind the observed relation. 


3.6 Luminosity Profile 

An alternative method to classify galaxies is to use the shape 
of the luminosity profile as a function of radius. Following 
Vincent & Ryden (2005), we used fracDeV to divide the RC 
in four subsets: galaxies with 0 <= fracDeV < 0.1 (n < 
1.2) are “ex” galaxies, galaxies with 0.1 <= fracDeV < 
0.5 (1.2 < n < 2.0) are “ex/de” galaxies. Conversely, 
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“de/ex“ galaxies have 0.5 <= fracDeV < 0.9 (2.0 < n < 

3.3) and finally “de” galaxies have fracDeV >= 0.9 (n > 

3.3) . 

In Fig. 10 we plot the e-R relation for the RC, 
splitting the galaxies according to their light profile. For 
flat ex galaxies the best-fit slope shows that there is no 
evidence of correlation between e and log R (top panel). 
For flat ex/de galaxies there is marginal evidence for 
a trend (the slope of the best-fit linear relation is 2 a 
above 0, second panel). For flat de/ex and de galaxies the 
significance increases to 5 a (third and bottom panels). 
The fraction of flat ex and ex/de galaxies increases with 
log R (numbers in parentheses). This trend, combined 
with the fact that ex and ex/de galaxies are on average a 
little flatter than de/ex and de galaxies, contributes to the 
e-R relation observed for the undivided RC (see Section 3.1). 
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log R [kpc] 


Figure 10. Same as Fig. 6, but here we split the sample according 
to the shape of the luminosity profile. For flat galaxies, the trend 
of increasing e with log R is recovered for de/ex and de galaxies 
(third and fourth panel). For flat ex/de galaxies there is only 
marginal evidence of a trend (second panel), while no significant 
trend is observed for flat ex galaxies (top panel). See Section 3.6 
for the definition of the labels. 
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4 DISCUSSION 

We have shown that in rich clusters e depends on radius. We 
find that galaxies closer to the cluster centre are on average 
rounder. The e-R relation persists when considering only 
galaxies with e ^ 0.4. In the next section we evaluate the 
effect of observational bias and establish that the trend is 
indeed genuine. We then evaluate the effect of the luminosity 
(Section 4.2) and of the T-E relation (Section 4.3). Finally 
we discuss two related works (Section 4.4) and proceed to 
review the possible explanations for the observed relation 
(Section 4.5). 


4.1 Observation bias 

The value of e is already corrected for PSF effects, because 
the model for the light profile is convolved with the relevant 
PSF. To evaluate any residual PSF effect we conduct two 
tests. We repeated our analysis using e m0 ments instead 
of tmodei, (as in e.g. KR05). e m oments was corrected for 
PSF effects according to the prescription of Bernstein & 
Jarvis (2002). Given that e mo dei is systematically higher 
than Cmoments , it is not surprising that our result changes. 
However the observed trends are still present (see Appendix 
A). As a further test for the effect of the PSF, we repeated 
our analysis limiting our sample to well resolved objects 
(where t/tpsf > 6.25, see KR05). The e-R relation persists 
even for these galaxies where the effect of the PSF is 
negligible, with lower significance as expected from the 
smaller sample size. 

Finally, we remark that the PSF would make all 
galaxies appear rounder, affecting smaller galaxies more 
than larger ones. This means that in order for PSF effects 
to cause or enhance the e-R relation, smaller galaxies would 
need to be a larger fraction of the population at small 
R, but the opposite is true (Dressier 1980; Hogg et al. 2004). 

All of our samples include a number of interlopers, fore¬ 
ground or background galaxies which are not gravitation¬ 
ally bound to the relevant cluster. This is true both for the 
parent sample and (to a lesser extent) for the RS sample. 
To understand the possible effect of interlopers on the ob¬ 
served trend of e versus log R, we consider the following 
scenario. For a given cluster, we assume that interlopers are 
distributed uniformly across the sky. In contrast, the den¬ 
sity of cluster member galaxies is higher in the centre and 
decreases at larger radii. This means that the fraction / of 
interlopers over the sample increases with log R. If interlop¬ 
ers had a higher average e than cluster members, then the 
observed variation of e with log R would be a consequence 
of the variation of / with log R. For Abell 1656 this sce¬ 
nario is ruled out by the fact that the e-R relation holds 
even for the sample of spectroscopically confirmed cluster 
members (Section 3.2). In addition, the 61 interlopers have 
an average ellipticity e = 0.34 ± 0.19, lower than the value 
for cluster members. As expected, the fraction of interlopers 
increases with log R. If interlopers have a lower than aver¬ 
age e for the other clusters too, then their inclusion cannot 
artificially create the observed e-R relation, but rather acts 
to mask it. 


4.2 The effect of luminosity 

More luminous galaxies are both rounder and more com¬ 
mon in the densest environments (Dressier 1980; Hogg et al. 
2004). Could the e-R relation arise from this trend? This 
hypothesis is contradicted by the fact that for the flat sub¬ 
sample there is no trend of e with M r (Fig. 8), but for the 
same galaxies a clear correlation exists between e and log R. 
In addition, the trend of e with log R is observed even at 
fixed magnitude, and is stronger for the less luminous galax¬ 
ies (Table 2). 


4.3 The morphology-density relations 

The final source of bias that we consider is the T-E relation 
(Dressier 1980). Es, SOs and spiral galaxies have different e 
distributions. Therefore the trend of their number fraction 
with local environment affects systematically the overall e 
distribution. In order to quantify this effect, we need Hubble 
morphological classifications for all the galaxies, but these 
are not available. 

In Section 3.3 we divided the RC (and the reference 
RS) into a flat (e Js 0.4) and a round subsample (e < 0.4). 
The subset of round galaxies contains both SRs (regardless 
of their inclination on the sky) and close-to-face-on discs. 
Therefore the observed trend of increasing flattening with 
log R can be due to the decreasing fraction of SRs (Cap- 
pellari et al. 2011b; D’Eugenio et al. 2013; Houghton et al. 
2013; Scott et al. 2014; Fogarty et al. 2014) and/or to a 
change in the intrinsic shape of galaxies. Without a dynam¬ 
ical classification for each galaxy it is impossible to draw 
any conclusions from the observed trend. 

However the subset of galaxies with e Js 0.4 does not 
contain SRs (which are intrinsically round). Therefore the 
observed trend of increasing projected ellipticity with radius 
arises because flat galaxies have lower intrinsic ellipticity 
nearer the centre of clusters. Flat galaxies include both FR 
ETGs and spiral galaxies. If these two classes have different 
n(e), the e-R relation could result from a trend of the mor¬ 
phological fractions with cluster-centric radius. However, it 
seems that on average FR ETGs have the same intrinsic 
flattening as spiral galaxies (Weijmans et al. 2014), a fact 
that undermines the above explanation. Nevertheless in Sec¬ 
tion 3.5 we attempted to remove spiral galaxies: we used 
GZ2 data to select a sample of smooth objects. When we 
adopt the condition P sm0 oth > 0.8, we find 252 galaxies, 
of which only 50 are flat (e ^ 0.4). The slope of e-R re¬ 
lation for flat galaxies is 1.8 standard deviations above 0 
(Fig. 9, top panel). For a sample of 50 flat galaxies selected 
randomly from the RC, we expect the slope to be 0.6 stan¬ 
dard deviations away from 0. This suggests (but does not 
demonstrate) that the low significance of the e-R relation 
of smooth, flat galaxies is due to insufficient sample size, 
as opposed to the systematic removal of the galaxies with 
^smooth ^ 0.2. We then relaxed the morphological selec¬ 
tion, finding 676 galaxies, of which 204 are flat (we used 
Ts m ooth > 0.5 and P sp irai < 0.5, see again Section 3.5). In 
this case the slope of the e-R relation of flat galaxies was 2.1 
standard deviations above 0. For a sample of 204 flat galax¬ 
ies randomly selected from the RC, we expected the slope to 
be 1.3 standard deviations above 0. The fact that we find a 
value of 2.1 suggests again that the low significance is due to 
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the sample size rather than to any systematic effect. In ad¬ 
dition, the slope is statistically consistent with what found 
in Fig. 6 , even though a significant number of spiral galaxies 
have been removed. 

Willett et al. ( 2013) cross-correlated the GZ2 catalogue 
with a professionally classified sample of SDSS galaxies 
(Nair & Abraham 2010). They showed that galaxies classi¬ 
fied as smooth in GZ2 (P smoo th > 0.8) are morphological 
ETGs for 96.7 % of the cases (the matching for our relaxed 
conditions is not covered). Therefore it appears that the 
spiral contamination is not the driving mechanism behind 
the observed relation. 

Our result could still be explained by the T-E relation, 
because within ETGs the fraction of Es increases with 
projected density (and nearer to the centre of clusters, 
Whitmore et al. 1993). Since Es are on average rounder 
than SOs, the e-R relation could be an effect of the T-E 
relation. To investigate this possibility we need to distin¬ 
guish between Es and SOs, which cannot be done with 
GZ2 and is beyond the scope of this paper. But the latest 
view on the E/SO divide, offered by the ATLAS 3D survey, 
is that the traditional visual E/SO morphology is flawed. 
Emsellem et al. (2011) demonstrated that two thirds of the 
locally classified Es are kinematically identical to SOs. We 
therefore reframe our findings within the FR/SR paradigm. 
15% of the local ETGs are SRs, but their fraction increases 
dramatically nearer the centre of clusters. (Cappellari et al. 
2011b; D’Eugenio et al. 2013; Houghton et al. 2013; Scott 
et al. 2014; Fogarty et al. 2014). Since SRs are on average 
rounder than FRs, the e-R relation could be a consequence 
of the kinematic morphology-density relation. However, 
by selecting galaxies with e 0.4, we eliminate all SRs. 
By further imposing P smoo th > 0.8 we eliminate all spiral 
galaxies. The fact that the e-R relation persists even for 
smooth galaxies flatter than e = 0.4 demonstrates that 
the e-R relation is not a consequence of the kinematic 
morphology-density relation. The decrease of projected 
ellipticity e nearer the centre of clusters implies that the 
distribution of intrinsic ellipticity changes with log R. 
This is evidence of the effect of the cluster and/or local 
environment on the shape of FRs, be they FR Es or FR SOs. 


4.4 Comparison with other work 

KR05 used a much larger sample of « 200,000 galaxies to 
investigate the relation between the shape of galaxies (as 
expressed by the axis ratio q) and the local environment 
(measured by the number density of neighbour galaxies in 
the sky). Their sample differs from ours in that: (i) they 
adopted a cut in apparent magnitude at r' = 17.77 mag and 
(ii) their selection is not limited to specific environments, 
whereas we focussed on rich clusters of galaxies. These dif¬ 
ferences make it difficult to compare our results. 

The noise introduced by their selection can be appreci¬ 
ated by looking at their Tables 3-5. Despite having « 20 x 
the sample size, the significance of the recovered trends 
is similar to ours. KR05 divide their sample according 
to the luminosity profile, as we did in Section 3.6. Our 
results are in qualitative agreement with theirs. However 


they do not investigate the trend of intrinsically flat galaxies. 

Weijmans et al. (2014) studied the intrinsic shape dis¬ 
tribution of the ATLAS 30 ETGs, using Integral Field Spec¬ 
troscopy to separate SRs and FRs. As we have seen, they 
find that the intrinsic flattening of FR ETGs is consistent 
with that of spiral galaxies, a fact that reinforces our re¬ 
sult. However, contrary to our findings, they do not observe 
any correlation of intrinsic flattening and environment (as 
parametrised by log E 3 ). It is unclear if this depends on the 
small size of their sample, and/or on the range of environ¬ 
ments that they explore (the ATLAS 30 survey contains just 
one unrelaxed cluster, Virgo). 

4.5 Possible driving mechanisms 

Here we propose three possible explanations for the observed 
trend. 

(a) We can view our result in the framework of the B/D de¬ 
composition. By definition discs are intrinsically flat, while 
bulges are rounder. A systematic change in the B/D mass 
ratio with log R could therefore account for the e-R relation. 
We notice however that in their study of eight nearby clus¬ 
ters (z < 0.06), Hudson et al. (2010) do not find any correla¬ 
tion between R-band B/D light ratios and the cluster-centric 
radius (except for the innermost « 0.3 Mpc). 

(b) The cluster environment affects the stellar population of 
member galaxies, for example by removing cold gas from the 
disc and halting star formation there (see Boselli & Gavazzi 
2006, for a review). We therefore expect to observe a trend 
between the cluster-centric radius and the photometric prop¬ 
erties of galaxies. For instance, even at fixed B/D mass ratio, 
a trend in stellar age can generate a trend in the B/D light 
ratio, which in turn could give rise to the e-R relation. This 
is supported by the observed correlation between disc colour 
and cluster centric radius, while no such correlation is ob¬ 
served for bulge colour (Balogh et al. 2000; Hudson et al. 
2010; Head et al. 2014). 

(c) FRs might contain subclasses with different intrinsic 
ellipticity, for instance Weijmans et al. (2014) suggests the 
existence of a tail of rounder FRs in the ATLAS 30 sample. 
If this is correct, the e-R relation for FRs could arise from 
the radial distribution of different classes of FRs. 

(d) Cluster environments also affect the dynamical struc¬ 
ture of galaxies. Harassment and stripping make discs 
smaller and thicker, therefore changing their intrinsic ellip¬ 
ticity. This is supported by Houghton et al. (2013), who 
studied a sample of galaxies in the central 15' (radius) of 
the Coma Cluster. For each galaxy they measure A (a proxy 
for the projected angular momentum per unit mass Em¬ 
sellem et al. 2007), and find that the maximum value of A 
for Coma FRs is lower than that of FRs from the ATLAS 30 
survey (Emsellem et al. 2011). They suggest that this is due 
to the cluster environment affecting the anisotropy (Binney 
& Tremaine 1988) of galaxies. 

In order to understand the origin of the e-R relation, 
an extended integral field spectroscopic survey of cluster 
galaxies is required. The Coma Cluster is an ideal start¬ 
ing point. Ongoing surveys (SAMI, MaNGA; Groom et al. 
2012; Bundy et al. 2015) could also address this problem, 
depending on their coverage of rich clusters. 
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5 CONCLUSION 

We studied the projected ellipticity t of galaxies in a sample 
of twenty nearby (z < 0.1) rich clusters and show that: 

(i) there is no evidence of differences between the distribu¬ 
tion of e among the clusters in our sample, except possibly 
for the Coma Cluster (Abell 1656) 

(ii) there exists a correlation between t and the projected 
cluster-centric radius R , which we refer to as the t-R relation 

(iii) the t-R relation exists independently of the trend of 
more luminous galaxies to be both rounder and more com¬ 
mon in the centre of clusters 

(iv) the t-R relation is stronger for galaxies on the RS, and 
persists for a redshift selected sample in the Coma Cluster 

(v) the t-R relation is steeper for flat galaxies (e Js 0.4; 
excludes SRs) than for round galaxies (e < 0.4; a mixture of 
FRs and SRs) 

(vi) for flat galaxies, there is no correlation between t and 
r'-band luminosity 

(vii) there is marginal evidence (2 a) that the t-R relation 
persists for flat ETCs, as classified by GZ2. The low signif¬ 
icance is likely due to the small sample size. 

We concluded that the t-R is evidence for physical effects 
causing intrinsically flat galaxies (including FR ETCs) to 
be rounder near the centre of rich clusters. 
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APPENDIX A: MOMENTS ELLIPTICITY 


KR05 use the weighted second order moments of the r'- 
band surface brightness I to measure the shape of galaxies. 
Following Bernstein & Jarvis (2002) they use a Gaussian 
weight function w(r, c) matched to the size and shape of the 
galaxy being measured (r and c are the rows and columns 
in the image). If we dehne: 


= f f(r, c)w(r, c)/(r, c)dr dc 
w f w(r, c)I(r, c)dr dc 

then the centre of the galaxy is given by: 

(r 0 ,c 0 ) = {{r) w , (c) w ) 

and the (central) second moments are defined as: 

Mrr = ((r - r 0 ) 2 ) w 
M rc = ((r - r 0 )(c - Co)) w 
M cc = ((c - co) 2 )™ 

From these quantities we can define: 

t = M cc + M rr —> mRrCr 
M cc - M r - 


(Al) 


(A2) 


(A3) 


e+ 


■ rl 


2 Mr, 


mE2 


(A4) 


where mRrCc, mEl and mE2 are the corresponding entries 
in the SDSS database. The axis ratio is given by: 


9 


f 

1 + e 


and finally we define the projected ellipticity e as: 

e = l — q 


(A5) 


(A6) 


In order to correct for the effect of the PSF, we follow 
again Bernstein & Jarvis (2002), which use the fourth order 
moment as well as the moments of the PSF itself (this data 
is included in the SDSS database). 



log R [kpc] 

Figure Al. Same as Fig. 6, but using the projected moments 
ellipticity e m0 ments- The method of moments, as implemented in 
SDSS, gives smaller values than e mo del- The sample of flat galaxies 
(e ^ 0.4) is therefore smaller, and it is impossible to assess any 
trend from the GZ2 cross-correlated sample. 


When we use e = e m oments the t-R relation is still ob¬ 
served (Fig. Al). However, given that generally we have 
Emoments < tmodei (Section 2.1), we have fewer flat galax¬ 
ies (emoments ^ 0.4). Therefore we cannot use Emoments to 
constrain the t-R relation of smooth, flat galaxies (bottom 
panel in Fig. Al). 

































